Coupled charge and valley excitations in graphene quantum Hall ferromagnets 
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Graphene is a two-dimensional carbon material with a honeycomb lattice and Dirac-type low- 
energy spectrum. In a strong magnetic field, where Coulomb interactions dominate against disorder 
broadening, quantum Hall ferromagnetic states realize at integer fillings. Extending the quantum 
Hall ferromagnetism to the fractional filling case of massless Dirac fermions, we study the elementally 
charge excitations which couple with the valley degrees of freedom (so-called valley skyrmions). 
With the use of the density matrix renomalization group (DMRG) method, the excitation gaps are 
calculated and extrapolated to the thermodynamic limit. These results exhibit numerical evidences 
and criterions of the skyrmion excitations in graphene. 
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I. INTRODUCTION 

A recent experimental realization of single-layer 
graphene sheets 1 - has made it possible to confirm a num- 
ber of theoretical predictions of intriguing electric proper- 
ties of massless Dirac fermion systems, 2 including uncon- 
ventional quantum Hall effects (QHE)&^ with the half- 
integer Hall conductivity^ 



Ae\ 1. 



(1) 



at v = ±2, ±6, ±10, • ■ • , where a factor 4 is the Lan- 
dau level (LL) degeneracy, accounting for spin and valley 
symmetry in graphene. v — 2irl 2 B p is the filling factor, 
Ib = y/h/eB is the magnetic length, p is the carrier 
density measured from charge neutral Dirac point. 

In addition to these unconventional quantized values, 
which can be solely understood on the basis of the mass- 
less Dirac fermion spectrum in a magnetic field: 



E n = sgn(n)hv F < 



(2) 



recent experimental studies in a sufficiently strong mag- 
netic field revealed new quantum Hall states at v = 
0,±1,±4)£ where the electron-electron interaction may 
play a crucial role. Here relevant energy scales in 
graphene in a magnetic field are (i) Landau level (LL) 
separation between n = and n = ±1, hiu = 
y/2hv F /i B ^ 400 y/B[T] [K], (ii) Zeeman coupling, A 2 = 
gUB |B| ~ 1.5 x {B[ T])[K] , and (hi) the Coulomb en- 
ergy, e 2 /el B ^ 100 y 'B[T][K]. The activation energy 
measurements^ have shown that at v = ±4 the gap has 
linear B dependence and reasonably corresponds to A z , 
indicating Zeeman spin splitting. At v = ±1, however, 
the gap is approximately scaled by V~B, which indicates 
that the gap originates from the Coulomb interaction. 
There have been a number of theoretical investigations 
of these states, and there are two leading theoretical sce- 
narios for the origin of the gap. One is the quantum 
Hall ferromagnetism (QHF) 7 in which valley degrees of 
freedom, referred to as pseudospins, spontaneously split 



via the exchange energy at all integer fillingsi^ ' 10 ! 11 ! 12 ! 13 
Second is the spontaneous mass generation ; 14 ' 15 ! 16 which 
predicts a gap only in the n = LL. Although they 
are not entirely orthogonal to each other and the order 
parameters in the two theories can coexist^ excitation 
properties, including the existence of gaps in n ^ LLs, 
are different. 

In the QHF theory, the ground state at v n = 1, where 
v n is the filling factor for 7ith Landau level defined by 



v n =v-A{n- 1/2), 



(3) 



is fully spin and valley polarized, and the wave function 
can be represented by 



l^" = 1 }=II C U0), 



(4) 



where we assign the valley K and K' in graphene as 
z-component of pseudospin r = K or K' . Excita- 
tions from the symmetry broken states are described 
by (pseudo)spin wave and (pseudo)spin textures called 
skyrmion a 18 i 19 ' 20 i 21 i 22 or other types, depending on the 
LL index n. Yang et al*2. estimated skyrmion excita- 
tion energy within the framework of Hartree-Fock (HF) 
approximation. 

The fractional quantum Hall effects (FQHE) 
in graphene have been studied by Apalkov and 
Chakraborty 23 . As a consequence of the relativistic 
nature of electrons in graphene, the effective electron- 
electron interactions in n ^ LLs differ from that of 
conventional two-dimensional systems i 8 ' 11 ! 23 With the 
use of the exact diagonalizations they have calculated 
the magnetoroton energy at v n ~ 1/3 when spin and 
valley are fully polarized. Their results show that the 
magnetoroton energy in the n = 1 LL is larger than 
that in the n = LL in finite systems^ Nevertheless, 
the magnetoroton excitation is charge neutral and it is 
not connected to the charged gap which is necessary 
to realize the FQHE. Moreover, the valley degrees of 
freedom which couple with charge excitations must be 
taken into account since there is no external symmetry 
breaking effect for valley degrees of freedom. The 
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FIG. 1: The pseudo-spin (valley) polarized excitation gap A c 
and the pseudospin (valley) unpolarized (skyrmion) excitation 
gap A 3 at v n — 1/3 in the n — and 1 LLs. 



excitation structures and the size of the activation gaps 
in the thermodynamic limit have not yet been clarified 
systematically.— 

In this paper we extend the quantum Hall ferromag- 
netism to the fractional filling case of graphene and show 
analogous properties of pseudospin ferromagnetism and 
topological excitations in the fractional quantum Hall ef- 
fect (FQHE) states in graphene. We start with the pro- 
jected Coulomb interaction Hamiltonian onto a certain 
LL, and study charge and valley excitations, where we 
treat the valley degrees of freedom K and K' in the lan- 
guage of the pseudospin, while real spin degrees of free- 
dom are supposed to be frozen by the Zeeman splitting. 
We calculate the exact wave functions of QHF ground 
states and low-energy excited states basing on the den- 
sity matrix renormalization group (DMRG) method, and 
confirm the existence of the skyrmion excitations of val- 
ley degrees of freedom not only at integer filling v = 1 
but also at fractional filling v = 1/3 for the LL index 
n = 0, 1 and 2. The skyrmion energies at these fractions 
are calculated and extrapolated to the thermodynamic 
limit. The existence of the finite gaps at v n = 1 in the 
n 7^ LLs strongly supports the QHF but not the mass 
generation scenarios. We note that the DMRG method 
is suitable for the present study, since we need to treat 
a large number of basis of the many-body Hilbert space 
including the psudospin degrees of freedom. To the best 
of our knowledge, there is no numerical work so far to 
estimate reliable extrapolated energies of the skyrmion 
excitation in the thermodynamic limit even in the con- 
ventional FQHE systems^ 

II. MODEL AND METHOD 

We assume that the total magnetic field is strong 
enough to cause complete spin splitting, and neglect spin 



flip excitations. We also neglect LL mixing for simplicity, 
although the ratio (e 2 /ds)/(?iw ) ~ 0.25 is not so small. 
The effect of LL mixing is discussed later. Taking ac- 
count of valley degrees of freedom, we apply the DMRG 
methocl 2 ^ 2 -^ on the spherical geometry. 29 ' 30 ^! 3 ^ 3 . 

The projected Hamiltonian onto the nth Landau level 
is written as 3 - - 

^ w =EE^h, (5) 

where Py [m] projects onto states in which particles i and 

j have relative angular momentum Tim, and Vm is their 
interaction energy in the nth Landau level^ Using the 
relativistic form factor in the nth Landau level 8 - 

F R (q) = L {q 2 /2) (6) 

and 

F£ («) = (1/2) [L H (<z 2 /2) + L W _ X (g 2 /2)] , (7) 
the pseudopotential o 29 ' 30 are given by 

Vi n) = r %V{q)e-<? [F^(q)} 2 L m (q 2 ). (8) 
Jo 27r 

Here L n (x) are the Laguerre polynomials. The corre- 
sponding integrals for electrons on the surface of a sphere 
which are used in the present work are described in 
Refs. [30|l3l|l3l. Note that our Hamiltonian Eq.([5]) has 
SU(2) symmetry in the valley degrees of freedom. A sym- 
metry breaking correction to Eq.© which stems from 
the honeycomb lattice structure of graphene 9 -* 1 ^ is order 
of a/is (a being a lattice spacing) in units of e 2 /c£b and 
neglected in the following. 

We calculate the ground state wave function using the 
DMRG method, 26 which is a real space renormalization 
group method combined with the exact diagonalization 
method. The DMRG method provides the low-energy 
eigenvalues and corresponding eigenvectors of the Hamil- 
tonian within a restricted number of basis states. The 
accuracy of the results is systematically controlled by 
the truncation error, which is smaller than 10 -4 in the 
present calculation. We investigate systems of various 
sizes with up to 40 electrons in the unit cell keeping 1400 
basis in each block i 27 ' 28 

In the sphere geometry, the pseudospin (valley) polar- 
ized ground state at v = 1/q (q being an odd integer), 
the Laughlin stated realizes when the total flux is 
given by 3 ^ 

N (f> (p,N e ) = v- 1 (N e -l), (9) 

where N e is the number of electrons in the system. Ele- 
mentary charged excitations from this pseudospin polar- 
ized ground state correspond to the ground state configu- 
rations of the system with additional/missing flux ±1. In 
the following, we study two types of excitations: Laugh- 
lin's quasiholes (quasiparticles)^ 9 - and skyrmion quasi- 
holes (quasiparticles)i 18 i 19 ' 20 i 21 ' 2 » 2 Laughlin's quasiholes 
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FIG. 2: The pseudospin unpolarized (skyrmion) excitation 
gap A s at v n — 1 in the n = 0, 1 and 2 LLs. The crosses 
on the vertical axis represent the results obtained by HF 
calculations.— The pseudospin polarized excitation gaps in 
n = 2 LL are plotted by open circles. 



(quasiparticles) correspond to the (pseudo)spin polarized 
excitations with ±1 flux, whose creation energy is given 
by 

Af = E{N 4> ±1,P = l)-E{N^P = 1), (10) 

where ± represents quasiholes and quasiparticles, respec- 
tively, and P is the polarization ratio of the pseudospin, 
i.e., P = {N K -N K >)/(N K +N K <) with N K [N K >) being 
the number of electrons in K {K') valley. 

Skyrmion quasiholes (quasiparticles) correspond to the 
(pseudo)spin singlet excitations, and their creation en- 
ergy is given by 

Af = E(Nt±l,P = 0)-E(Nt,P=l), (11) 

which could be smaller than A^. 

The activation energy, referred to as the gap in the 
following, is given as a sum of these quasihole and quasi- 
particle energies, A c = A^+A^ for pseudospin polarized 
excitations, and A s = A+ + Aj for pseudospin unpolar- 
ized excitations. 



III. RESULTS 

Figure 1 shows A c and A s at v n = 1/3 as a function of 
1/N e . In the n = LL, the pseudospin (valley) polarized 
excitation gap A c is 0.101 e 2 /(e/s) in the thermodynamic 
limit in good agreement with the previous workj^L In the 
n = 1 LL, A c is 0.115, which is larger than 0.101 in 
the n = LL. This enhancement of A c stems from the 
unique properties of the pseudopotentials in n = 1 LL for 
relativistic particles; > V^ 0) , 23 The magneto- roton 
excitation energy in the n = 1 LL is also larger than that 
of the n = LL^ 



FIG. 3: Two-particle correlation functions g TT i of quasihole 
skyrmion state at v\ — 1/3. 



The larger excitation gap in the higher LL has also 
been obtained for the pseudospin (valley) unpolarized 
excitations A s ; the unpolarized excitation gap A s in the 
n = 1 LL is 0.05 e 2 / (els), which is larger than that in the 
n = LL, A s = 0.03. We note that similar analysis has 
been done in Ref. with qualitatively different results. 

The results for v n = 1 are shown in Fig. 2. The ex- 
trapolated value of the pseudospin unpolarized excita- 
tion gap A s for n = LL is 0.63 e 2 /(ds), which is 
consistent with the previous HF studies and exact diag- 
onalization studies^ Even in the n = 1 and n = 2 LLs, 
A s are almost the same with the HF results obtained in 
the four-component model for grapheneJ^ These results 
show that the HF trial states for skyrmion excitation is 
essentially correct. 

In the above analysis, we have neglected the effect of 
LL mixing. When we study higher LLs, we need to take 
account of LL mixing, because LL spacing of graphene 
decreases with the increase in its index n. In conven- 
tional two-dimensional systems, it has been argued that 
LL mixing reduces the activation gap in the fractional 
QHE regimes^ 

To study the skyrmion-like structure in the pseu- 
dospin unpolarized excited states, we have calculated the 
two-particle correlation functions ff TT /(r) i 22 ' 29 i 30 Figure 3 
shows g T T'{f) for the pseudospin unpolarized quasihole 
state at v\ = 1/3. We And clear peak structure around 
the origin r = for the electrons in the same valley, 
9kk{t)- On the other hand, the maximum appears at 
the opposite side on the sphere (r = ttR) for the elec- 
trons in different valleys, gKK'(r)- These structures are 
consistent with the skyrmion-like pseudospin structure. 

The HF trial state of quasihole skyrmions at v n = 1 is 
written in the form: 

l*«fc)= II [ a m4nK+0 m 4n+lK>¥), (12) 
m=-N 4 ,/2 
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FIG. 4: Expectation values, {cJ„ T c mT ), in pseudospin (valley) 
unpolarized excited states with one extra flux at v n =1 and 
1/3. 
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FIG. 5: The polarization ratio P dependence of the quasihole 
skyrmion energy AE(P) = E(Nt + 1, P) - E{N^ + 1, P = 0) 
at v n = 1 in the n = 0, 1 and 2 LLs. 



where <^ mK c mK ) = \a m \ 2 and {c X mK ,c mK ) = |/3 m _i| 2 . 
To see the relation to this trial state, we have calculated 
the expectation values (cJ rjr c mT ) from the wave func- 
tion obtained in the present DMRG study. The results 
shown in Fig. 4 indicate they are approximately given by 
(cJ nr c mT ) = 1/2 =p m/Ntj) for t = K and r = if', respec- 
tively. Similar results are also obtained for v n = 1/3 in 
the n = and 1 LLs. 

The skyrmion excitations at v n = 1 arc stable up 
to n = 2 LLs, as shown in Fig. 5, in which the polar- 
ization energies of the quasihole skyrmions AE(P) = 
^(iV + 1, P) - E(N<j> + 1, P = 0) are plotted. The clear 



P 2 -dependence shown in Fig. 5 and its inset confirms that 
the pseudospin singlet state is the lowest excitation for 
n = 0, 1 and 2 LLs. In the n — 2 LL, the minimum at 
P = appears only for systems whose number of elec- 
trons is larger than 16; that means the skyrmion excita- 
tions are stable only for large systems. This is also shown 
in Fig. 2 as an upward cusp in A s at N e ~ 18; the pseu- 
dospin singlet state becomes the lowest when the number 
of electrons exceeds 18, while the pseudospin polarized 
states marked by the open circles become the lowest for 
systems with N e < 18. 

In higher LLs (n > 3), pseudospin singlet excitations 
have not been seen as the lowest excitation within our 
calculations up to 40 electrons. These results are quali- 
tatively consistent with the HF analysis in Ref. 10, which 
predicts the skyrmion excitation gap is smaller than the 
HF quasiparticle gap only for lower LLs below n = 3. 
The skyrmion excitation gap in n = 3 LL is close to 
the HF quasiparticle gap and it is difficult to stabilize 
skyrmions in n = 3 LL. 



IV. DISCUSSION 

Our DMRG calculation confirms the valley polarized 
ground state at v n = 1 for n < 2 LLs and at v n = 1 /3 for 
n < 1 LLs. The (pseudo)spin polarized ground state at 
v n = 1/q can be written simply as Laughlin's Jastrow 
function^ and elementally charge excitations are ob- 
tained by increasing or decreasing the flux quantum num- 
ber N,p by 1. We have studied both (a) pseudospin po- 
larized excitations (Laughlin's quasiholes and quasiparti- 
cles) and (b) pseudospin unpolarized excitations (quasi- 
hole skyrmions and quasiparticle skyrmions). The acti- 
vation energies obtained in finite systems are extrapo- 
lated to the thermodynamic limit, which give theoretical 
predictions for future experimental studies of the frac- 
tional quantum Hall states in graphene. In a high qual- 
ity graphene sample, the v n = 1 QH states have been 
observed in the n = LL;S Very recently, suspended 
graphene sheets with metallic contacts were succeeded^ 
The reported mobility 2 x 10 5 cm 2 /Vs is close to the 
crudely obtained critical value where the FQHE real- 
izes within available magnetic field_< 50 T, estimated 
from the theory described in Ref. [8| and the skyrmion 
gap A s ~ 0.05(e 2 /els) obtained in the present work. 
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